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2 ■ Abstract 

A model describing Ising spins with short range interactions moving ran- 
fl_j ' domly in a plane is considered. In the presence of a hard core repulsion, 



which prevents the Ising spins from overlapping, the model is analogous to a 



p\ ' dynamically triangulated Ising model with spins constrained to move on a flat 



surface. As a function of coupling strength and hard core repulsion the model 
exhibits multicritical behavior, with first and second order transition lines 
terminating at a tricritical point. The thermal and magnetic exponents com- 
puted at the tricritical point are consistent with the KPZ values associated 
with Ising spins, and with the exact two-matrix model solution of the ran- 
dom Ising model, introduced previously to describe the effects of fluctuating 
geometries. 



I. INTRODUCTION 

Following the exact solution of the Ising model on a random surface by matrix model 
methods |]l[ , there has been growing interest in the properties of random Ising spins coupled 
to two-dimensional gravity. More recently, work based on both series expansions |^,^ and 
numerical simulations Q,^] has verified and extended the original results. It is characteristic 
of these Ising models that the spins are allowed to move at random on a discretized version 
of a fluid surface. In a specific implementation of the model, Ising spins are placed at the 
vertices of a lattice built out of equilateral triangles, and the lattice geometry is then allowed 
to fluctuate by varying the local coordination number through a "link flip" operation which 
varies the local connectivity [Q. Remarkably, the same critical exponents have also been 
found using consistency conditions derived from conformal field theory for central charge 
^ ~ 2 ' which should again apply to Ising spins. It is generally believed that the new values 
for the Ising critical exponents are due to the random fluctuations of the surface in which 
the spins are embedded, and therefore intimately tied to the intrinsic fractal properties of 
fluctuating geometries. It came therefore as a surprise that non-random Ising spins, placed 
on a randomly fluctuating geometry but with fixed spin coordination number, exhibited the 
same critical behavior as in flat space, without any observed "gravitational" shift of the 
exponents ^. 

The natural question is then to what extent the values of the critical Ising exponents 
found by KPZ for c = | and in the matrix model solution (a = —1, (3 = 1/2, ■y = 2, r] = 2/3, 
u = 3/2 ||l]J^) are due to the annea/ed randomness of the lattice, and to what extent they are 
due to the physical presence of a fluctuating background metric. The most straightforward 
way to answer this question is to investigate the critical properties of annealed random Ising 
spins, with interactions designed to mimic as closely as possible the dynamical triangulation 
model, but placed in flat two-dimensional space. It is well known that for a quenched 
random lattice the critical exponents are the same as on a regular lattice [§] , as expected on 
the basis of universality, even though in two dimensions the Harris criterion (which applies to 



quenched impurities only) does not give a clear prediction, since the specific heat exponent 
vanishes, a = 0, for Onsager's solution. 

In this paper we present some detailed results concerning the exponents of such a model 
in order to complete the discussion presented in a recent publication 0. 

II. FORMULATION OF THE MODEL 

In a square (i-dimensional box of sides L with periodic boundary conditions we introduce 
a set oi N = L'^ Ising spins Si = ±1 with coordinates x^, i = 1...N, a = l...d, and average 
density p = N/L'^ = 1. Both the spins and the coordinates will be considered as dynamical 
variables in this model. Interactions between the spins are determined by 

i[x, 5] = - E M^^^ ^j) w^j s^Sj -hJ2w,s, , (1) 

i<j i 

with ferromagnetic coupling 

{0 if \xi — Xj\ > R 
' \ ) 

J ii r < \xi — Xj\ < R 

and infinite energy for \xi — Xj\ < r, giving therefore a hard core repulsion radius equal to 
r/2. As will be discussed further below, the hard core repulsive interaction is necessary for 
obtaining a non-trivial phase diagram, and mimics the interaction found in the dynamical 
triangulation model, where the minimum distance between any two spins is restricted to be 
one lattice spacing. For r — > 0, Jij = J[l — 9{\xi — Xj\ — R)]. 

The weights Wij and Wi appearing in Eq. (|l|) could in principle contain geometric factors 
associated with the random lattice subtended by the points, and involve quantities such as 
the areas of the triangles associated with the vertices, as well as the lengths of the edges 
connecting the sites. In the following we will consider only the simplest case of unit weights, 
Wij = Wi = 1. On the basis of universality of critical behavior one would expect that 
the results should not be too sensitive to such a specific choice, which only alters the short 
distance details of the model, and should not affect the exponents. 



The full partition function for coordinates and spins is then written as 

N d L 

Z = I[ E (n/ dxt) exp(-/[x,5]) . (3) 

i=l 5,=±1 a=l -^0 

In the following we will only consider the two-dimensional case, d = 2, for which specific 
predictions are available from the work of KPZ and the matrix model solution. 

It should be clear that if the interaction range R is of order one, then, for sufficiently large 
hard core repulsion, r -^ v5/2 < R, the spins will tend to lock in into an almost regular 
triangular lattice. As will be shown below, in practice this crossover happens already for 
quite small values of r. The critical behavior is then the one expected for the regular 
Ising model in two dimensions, namely a continuous second order phase transition with the 
Onsager exponents. Indeed for the Ising model on a triangular lattice it is known that 
Jc = |v31n3 = 0.9514.... On the other hand if the hard core repulsion is very small, then 
for sufficiently low temperatures the spins will tend to form tight ordered clusters, in which 
each spin interacts with a large number of neighbors. As will be shown below, this clustering 
transition is rather sudden and strongly first order. Furthermore, where the two transition 
lines meet inside the phase diagram one would expect to find a tricritical point. 

In order to investigate these issues further, we have chosen to study the above system 
by numerical simulation, with both the spins and the coordinates updated by a standard 
Monte Carlo method. The computation of thermodynamic averages is quite time consuming 
in this model, since any spin can in principle interact with any other spin as long as they get 
sufficiently close together. As a consequence, a sweep through the lattice requires a number 
of order A^^ operations, which makes it increasingly difficult to study larger and larger 
lattices. In order to extend our study to even larger lattices, we have applied a binning 
procedure in such a way that the time for the updating of a given configuration grows as 
zN, where z is the average coordination number of the lattice, instead of A^^. This binning 
procedure consists of dividing the system in cells of unit length, and keeping track of the 
spins in each cell. Since all the moves are local, and spins can only move from a given cell to 
the neighboring ones, we only need to consider the spins in a given cell and its neighbors at 



each updating step. This procedure is very effective when the average coordination number 
is relatively small (J < Jc and r large), however, ii z ^ N the updating time grows again 

as Ar2. 

In this paper we will not address in detail the problem of critical slowing down, however 
an additional possibility for the future could be to implement some sort of cluster updating 
algorithm ]ll|. On the other hand, we should add that we have not found any anomalous 



behavior as far as the autocorrelation times are concerned, which remain quite comparable 
to the pure Ising case. 

There are a number of local averages and fluctuations which can be determined and 
used to compute the critical exponents. In the course of the simulation the spontaneous 
magnetization per spin 

was measured (here the averages involve both the x and S variables, < >=< >x,s'), as well 
as the zero field susceptibility 

^ = l^W-^^^\^-^ = l^<T.s-^s,>-]^<\Y.s^\>' . (5) 

It is customary to use the absolute value on the r.h.s., since on a finite lattice the spontaneous 
magnetization, defined without the absolute value, vanishes identically even at low temper- 
atures. In addition, in order to determine the latent heat and the specific heat exponent, 
we have computed the average Ising energy per spin defined here as 

Id 1 

and its fluctuation, 

Some additional quantities we have used in the course of this work will be defined later. 



III. RESULTS AND ANALYSIS 

In the simulations we have investigated lattice sizes varying from 5^ = 25 sites to 30^ = 
900 sites. The length of our runs varies in the critical region {J r^ J^) between 2M sweeps on 
the smaller lattices and 200k sweeps on the largest lattices. A standard binning procedure 
then leads to the errors reported in the figures. 

As it stands, the model contains three coupling parameters, the ferromagnetic coupling 
J, the interaction range R and the hard core repulsion parameter r. We have fixed R = 1; 
comparable choices should not change the universality class. As we alluded previously, 
for small r we find that the system undergoes a sharp first order transition, between the 
disordered phase and a phase in which all spins form a few very tight magnetized clusters, in 
which the number of neighbors is of the order N. These clusters persist even for larger values 
of the hard core repulsion, r, but the number of interacting neighbors does not become as 
large as A^ in this case. 

In Figs, in and |^ we show the existence of these clusters when the hard core repulsion 
is as large as r = 0.4. In Fig. |I| we observe ferromagnetic order in small domains even 
though we are below Jc- On the other hand, in Fig. ^ where we are above Jc, the system 
has practically clustered into a single ferromagnetic domain. For sufficiently large r, the 
transition is Ising-like, between ordered and disordered, almost regular, Ising lattices (for 
our choice of range R, the transition appears to be very close to regular Ising-like for r ^ 0.6 
and larger, see below). In Fig. ^ we show a particular configuration for r = 0.98 where the 
regular, almost triangular, lattice is clearly visible. In this case the average coordination 
number z is very close to 3, as expected for a regular triangular lattice. In Fig. ^ we show 
the average number of neighbors z for several values of r on a system with N = 144 spins. 
We find that for small values of r the coordination number increases very rapidly as we 
approach the critical point. On the other hand, for intermediate choices of r, z saturates to 
a smaller value. When r = 0.6, the coordination number saturates to a value oi z = 3.1, 
which is already very close to the value on a regular triangular lattice {z = 3). 



In Fig. ^ we plot the average energy per bond E^ as a function of J for several choices 
of the hard core repulsion r. The jump discontinuity, which is visible for small hard core 
repulsion r, indicates the existence of a first order transition. For larger values of r, the 
discontinuity is reduced and eventually vanishes. A determination of the discontinuity in 
the average energy of Fig. ^ at the critical coupling Jc shows that it gradually decreases as 
r is increased from zero. Fig. ^ shows a plot of the latent heat per bond A^ versus r at the 
transition point Jc- In general we do not expect the latent heat to vanish linearly at the 
endpoint, but our results seem to indicate a behavior quite close to linear. From the data 
we estimate that the latent heat vanishes at r = 0.344(7), thus signaling the presence of a 
tricritical point at the end of the first order transition line. Beyond this point, the transition 
stays second order, as will be discussed further below. The phase transition line is shown in 
Fig. |15|; for r = we found on the largest lattices Jc = 0.19(2), while for r = 0.98 we found 
Jc = 0.93(3). 

In Fig. 1^ we plot the spin susceptibility as a function of J for several system sizes near the 
tricritical point, showing a growth of the peak with system size. To determine the critical 
exponents, we will resort to a finite-size scaling analysis. In the following we will be mostly 
concerned with the values for the critical exponents in the vicinity of the tricritical point. 
In the case of the spin susceptibility, from finite-size scaling, we expect a scaling form of the 
type 

To recover the correct infinite volume result one needs x{^) ~ ^'"^ fo^' large arguments. 
Thus, in particular the peak in x should scale like N^/'^'^ for sufficiently large A^. In Fig. ^ 
we show the evolution of the computed peaks in x as a function of In A^. 

Despite the fact that the lattices are quite small, as can be seen from the graph, a linear 
fit to the data at the tricritical point is rather good, with relatively small deviations from 
linearity, x^ /d.o.f. ~ 10"^. Using least-squares one can estimate 7/z/. For r = 0.35 we find 
7/1/ = 1.32(3), which is much smaller than the exact regular Ising result 7/// = 1.75. From 



scaling one then obtains the anomalous dimension exponent rj = 2 — ^jv = 0.68(3). To 
further gauge our errors, we have computed the same exponent for the regular Ising limit, 
for r = 0.6. In this case we indeed recover the Onsager value: we find on the same size 
lattices and using the same analysis method 7/z/ = 1.72(4). We also note that the shift in 
the critical point on a finite lattice is expected to be determined by the correlation length 
exponent u, namely JdN) — Jc(oo) ~ N^^l"^^ . This relationship can be used to estimate z/, 
but it is not very accurate. From a fit to the known values of Jc{N) we obtain a first rough 
estimate z/ = 1.3(2). A more precise determination of u will be given later. 

A similar finite-size scaling analysis can be performed for the magnetization, which is 
shown in Fig. |^ for several system sizes. Close to and above Jc we expect M ~ {J — Jc)^- At 
the critical point on a finite lattice, as determined from the peak in the susceptibility (which 
incidentally is very close to the infiection point in the magnetization versus J), M should 
scale to zero as Mn(Jc) ~ A^^/^^. in Fig. [1^ we show the magnetization M computed in 
this way for different size lattices At the tricritical point we find jS/u = 0.31(4), which again 
clearly excludes the pure Ising exponent, /5/z/ = 0.125. For the pure Ising limit (r = 0.6) we 
obtain P/u = 0.13(7), which is close to the expected Onsager value. 

The results for the Ising specific heat C at the tricritical point as a function of lattice size 
A^ are shown in Fig. |ll]. One expects the peak to grow as C ~ iV"/^'", but the absence of any 



growth for the larger values of A^ implies that a/u < (a weak cusp in the specific heat). 
In general close to a critical point, the free energy can be decomposed into a regular and a 
singular part. In our case the singular part does not seem to be singular enough to emerge 
above the regular background, leading to an intrinsic uncertainty in the determination of an 
a < 0, and which can only be overcome by determining still higher derivatives of the free 
energy with respect to the coupling J. In order to isolate the singular part of the specific 
heat we have therefore calculated dC/dJ from the expression 

dC .„ 



17 = '' 



3{E){E')-{E')-2{Ey\. (9) 

In the infinite system dC/dJ should diverge according to 



§Hi-^.r'-'. (10) 



In particular, if a = —1, dC/dJ should diverge logarithmically. In Fig. |T^ we show the 
scaling of dC/dJ on a lattice with A^ = 256 spins according to Eq. (|10D- From the slope 
of the curve we determine the critical exponent to be a ~ —0.98(4). We have also tried 
to assume a logarithmic scaling behavior as shown in Fig. |1^. It is clear that from the 
linear behavior of dC/dJ we can conclude that our results are completely consistent with 



an exponent of a = —1. We attribute the small discrepancy between the results of Figs. ^ 
and [1^ to the fact that we are not sufficiently close to Jc and that we are on a finite lattice 
with A^ sites. We have also performed a similar analysis for the fiuctuation in the energy 
per bond (as opposed to the energy per site as defined previously). In this case we find close 
to the tricritical point a ~ —0.96(2). 

In the regular Ising case one has in a finite volume a logarithmic divergence C ~ In A^ 
(and q;/2z/ = 0), and we indeed see such a divergence clearly for r = 0.6, which corresponds 
to the almost regular triangular Ising case. 

Another approach to obtaining a is to determine the correlation length exponent u 
directly instead, and use scaling to relate it to a = 2 — 2z/. The exponent u can be obtained 
in the following way. First one can improve on the estimate for Jc by considering the 



fourth-order cumulant [|I2 



f/^(j) = 1 - -^^4^ , (11) 

where rn = J2i Si/N. It has the scaling form expected of a dimensionless quantity 

Un{J) = U{N'/'^\J-Jc\). (12) 

The curves 11^ (J), for different and sufficiently large values of A^, should then intersect at a 
common point J^, where the theory exhibits scale invariance, and U takes on the fixed point 



value U*. In Fig. 14 we show the fourth-order cumulant as a function of J for r = 0.35 and 



for several lattice sizes. We have found that indeed the curves meet very close to a common 

9 



point, and from the intersection of the curves for A^ = 25 to 400 we estimate Jc = 0.472(9), 
which is consistent with the estimate of the critical point derived from the location of the 
peak in the magnetic susceptibility. We also determine U* = 0.47(4), which should be 
compared to the pure Ising model estimate for the invariant charge U* ~ 0.613 [|T3|. 



One can then estimate the correlation length exponent z/ from the scaling of the slope of 
the cumulant at Jc- For two lattice sizes N, N' one computes the estimator 



with U'pf = dllN/dJ defined by 



rt, ^ 



'N 



3{m^y 



{m^){E) + {m^E)-2- 



rrr){m^E) 



(14) 



Using values of A^ from systems with 256, 400, and 900 spins we estimate v from Eq. ([M]) 
to be 1.46(8). Using the scaling relationship a = 2 — 2z/, we obtain an estimate for a which 
is again quite consistent with our previous estimate derived from dC/dJ. 

In Table |I| we summarize our results, together with the exponents obtained for the two- 
matrix model |I| (and which are the same as the KPZ values |0), for the Onsager solution 
of the square lattice Ising model, and for the tricritical Ising model in two dimensions [ITU 



As can be seen, the exponents are quite close to the matrix model values (the pure Ising 
exponents seem to be excluded by several standard deviations). 
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TABLES 
TABLE L Estimates of the critical exponents for the random two-dimensional Ising model, as 

obtained from finite-size scaling at the tricritical point. 





-f/v 


13/v 


a/u 


a 


V 


This work 


L32(3) 


0.31(4) 


-0.65(4) 


-0.98(4) 


1.46(8) 


Matrix model 


1.333... 


0.333... 


-0.666... 


-1.0 


1.5 


Onsager 


1.75 


0.125 


0.0 


0.0 


1.0 


Tricritical Ising 


1.85 


0.075 


1.60 


0.888... 


0.555... 
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IV. CONCLUSIONS 

In the previous sections we have presented some results for the exponents of a random 
Ising model in flat two-dimensional space. The model reproduces some of the features of a 
model for dynamically triangulated Ising spins, and in particular its random nature, but does 
not incorporate any effects due to curvature. Due to the non-local nature of the interactions 
of the spins, only relatively small systems have been considered so far, which is reflected in 
the still rather large uncertainties associated with the exponents. Still a rich phase diagram 
has emerged, with a tricritical point separating first from second order transition lines. The 



phase diagram we obtain is shown in Fig. |15[ We have localized the tricritical point at 



Jc = 0.471(5) and r = 0.344(7). The thermal and magnetic exponents determined in the 
vicinity of the tricritical point (presented in Table |ID have been found to be consistent, within 
errors, with the matrix model solution of the random Ising model and the KPZ values. Our 
results would therefore suggest that matrix model solutions can also be used to describe a 
class of annealed random systems in fiat space. 

One might wonder at this point if the spin system discussed in this paper can be found 
among the models in the FQS classification scheme [|1^] of two-dimensional conformally 
invariant field theories 0. Since the model is apparently not unitary (it contains short range 
repulsion and long range attraction terms), it should fall into the wider class of degenerate 
theories considered by BPZ |]15|. The allowed scaling dimensions in these theories are given 



by the well-known Kac formula, 

_ 1 

P'9 ~ 4 



(pa+ + qa^Y — (a+ + a_)^ (15) 



with p,q positive integers, and a± = ao ± (1 + a'^Y^'^. Oq is related to the conformal 
anomaly c of the theory by c = 1 — 24aQ. Often the central charge is then written as 
c = 1 — 6/m{m + 1). One of the difficulties in this approach is the identification of a 



^We thank Giorgio Parisi for suggesting to look into this aspect. 
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given realization of confornial symmetry with a particular universality class. The simplest 
possibility appears to be m = 4/5, corresponding to m = r/{s — r) with s = 9 and r = 4. 
One then obtains for this choice the central charge c = —19/6, and ao = 5/12, a+ = 3/2 
and a_ = —2/3. The matching scaling dimensions are then Ai^4 = As^s = 1/6 (which gives 
1] = 2/3), and Ai 5 = A3 4 = 2/3 (which gives u = 3/2). Negative values of c are allowed in 
non-unitary theories. It would be of interest to compute the central charge directly in the 
random spin model and verify this assignment, using the methods described in Ref. |T7|| . 

We should mention in closing that the above values for s and r appear to be rather close 
to the ones associated with the Yang-Lee edge singularity, which describes the behavior of 
the magnetization in the Ising model in the presence of an imaginary external field, and for 
which Cardy [1^ has suggested the identification s = 5 and r = 2, which yields m = 2/3 



and c = —22/5. It is known that the Yang-Lee edge singularity also describes the critical 
properties of large branched dilute polymers and of the Ising model in a quenched random 
external field in d + 2 dimension [|T8| . 
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FIGURES 
FIG. 1. A particular configuration of the system with N = 400 spins and r = 0.4 for J = 0.35. 

Spins with S = ±1 are indicated with empty and solid circles, respectively. 

FIG. 2. A particular configuration of the system with N = 400 spins and r = 0.4 for J = 0.65. 
Spins with S = ±1 are indicated with empty and solid circles, respectively. 

FIG. 3. A particular configuration of the system with N = 400 spins and r = 0.98 for 
J = 0.25. The hard core repulsion radius is shown as a circle around the spin. Spins with S = ±1 
are indicated with empty and solid circles, respectively. 

FIG. 4. The average number of neighbors z as a function of J on a lattice with N = 144 spins 
for several choices of the hard core repulsion r. 

FIG. 5. The average energy per bond Ez as a function of J for several choices of the hard core 
repulsion r for a system with N = 100 sites. 

FIG. 6. The latent heat per bond A^ along the first order transition line, plotted against the 
hard core repulsion parameter r. The tricritical point is located where the latent heat vanishes. 

FIG. 7. The magnetic susceptibility x versus J for fixed hard core repulsion parameter r = 0.35 
and different system sizes. 

FIG. 8. The peak in the magnetic susceptibility Xmax versus the number of Ising spins N for 
choices of the hard core repulsion parameter corresponding to r = 0.35 and r = 0.6. 

FIG. 9. The magnetization M versus J, for fixed hard core repulsion parameter r = 0.35 and 
different system sizes. The solid line is a spline through the data for N = 144. 

FIG. 10. Finite size scaling of the magnetization at the inflection point Mi^f versus the total 
number of Ising spins N for choices of the hard core repulsion parameter corresponding to r = 0.35 
and r = 0.6. 
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FIG. 11. Plot of the specific heat C versus ferromagnetic couphng J at r=0.35, showing the 
absence of a growth in the peak with increasing lattice size (for the larger systems) , in contrast to 
the behavior of the magnetic susceptibility. The errors (not shown) are smaller than the size of 
the symbols. 

FIG. 12. The derivative of the specific heat dC/dJ as a function of Jc — J on logarithmic axes 
for iV = 256. 

FIG. 13. The derivative of the specific heat dC/dJ as a function ol Jc — J on semi-logarithmic 
axes for A^ = 256. 

FIG. 14. The Binder fourth-order cumulant [/ as a function of J for fixed hard-core repulsion 
r = 0.35 and on several lattices with A^ spins. The solid line is a spline through the data for 

N = 144. 

FIG. 15. The phase diagram for the dynamical random Ising model on a two-dimensional flat 
surface. The tricritical point (denoted by the solid circle) separates the first order from the second 
order transition lines. The paramagnetic (PM) and ferromagnetic (FM) phases are also shown. 
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